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Abstract 


The experimental charge and discharge profiles of a LiCoO, electrode show that the overpotential of the electrode does not change much during 
galvanostatic charge, but changes significantly during galvanostatic discharge. Semi-empirical porous electrode models are presented to simulate 
the charge and discharge profiles of the LiCoO> electrode. The symmetry factor is empirically assumed to decrease with the state of discharge of 
the electrode to enable the model predictions to agree well with the experimental discharge profiles. 
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1. Introduction 


Lithium cobalt oxide has been widely used as cathode mate- 
rial in lithium ion secondary batteries. The reaction of LiCoO}? is 
well known to be lithium ion extraction from and insertion into 
a layered cobalt dioxide matrix. Many researchers [1-3] have 
simulated the charge and discharge behaviors of a full cell sand- 
wich which uses LiCoO2 as cathode material. However, very 
few studies can be found in the literature focusing on LiCoO2 
electrode alone. In this work, the rate capability of a LiCoO2 
electrode was studied in half cells with a three-electrode setup. 
The study shows that the overpotential of the LiCoO> electrode 
does not change much during galvanostatic charge, but changes 
significantly with cell state of discharge during galvanostatic 
discharge. A pseudo-2D porous electrode model [4—6] was used 
to simulate charge and discharge profiles of the LiCoO2 elec- 
trode. Empirical expressions were used in the Butler-Volmer 
equation to enable the model predictions to agree well with the 
experimental data. 
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2. Experiment 


The rate capability of a LiCoO2 electrode was measured using 
a Swagelok-type half cell (see Fig. 1). The half cell consisted 
of a LiCoOz working electrode, a separator and a lithium foil 
counter electrode. Cone shaped lithium metal with a sharp end 
was placed just above the LiCoO electrode to serve as the ref- 
erence electrode. The potential of the LiCoO2 electrode versus 
Li/Li* was measured directly in the experiment. The working 
electrode was a round disc with a diameter of 1/2 in. pouched 
out of a singled-sided LiCoOz2 sheet electrode provided by Mine 
Safety Appliances (Sparks, MD). The Celgard-2300 polypropy- 
lene membrane (Charlotte, NC) of a thickness of 25 ym was used 
as the separator. The electrolyte (Ferro, Independence, OH) was 
1.0M LiPF¢ in a solvent mixture of EC/PC/EMC/DEC. The 
half cell was assembled in an argon-filled glovebox, removed 
and placed in a Tenney environmental chamber controlled at 
15°C. An 8-channel Arbin battery test unit was used to conduct 
the rate capability test. 

The half cell was first cycled several times at the C/13 rate 
(C=4mA) between 3.0 and 4.35 V versus Li/Li* to stabilize 
the electrochemical performance of the LiCoOz electrode. The 
capacity of the LiCoO, disc electrode was measured around 
4mAha when it was cycled between the given voltage win- 
dow. The currents used in rate capability test included 6 x 1075, 
3 x 1074, 6 x 1074, 1.5 x 1073 and 3 x 107? A, which roughly 
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Nomenclature 
ds specific interfacial area of the electrode 
(cm? cm7?) 
Brug Bruggeman coefficient 
Ce Li* concentration in the electrolyte (mol cm~?) 
Cs solid phase Li* concentration in the LiCoO3 elec- 


trode (mol cm7?) 


Csmax maximum solid phase Li* concentration in the 
LiCoO electrode (mol cm~?) 

De diffusion coefficient of the electrolyte (cm? s7!) 

D; solid phase diffusion coefficient (cm? s7!) 

fa mean molar salt activity coefficient 

F Faraday’s constant, 96487 C mol™ ! 

io exchange current density (A cm7?) 

I applied current (A) 

j intercalation current density (A cm~7) 

kp kinetic rate constant (A cm7? (mol cm~+)~!5) 

Q charge or discharge capacity (Ah) 

R gas constant, 8.3145 J (mol kK)! 

Rf contact resistance (Q cm”) 

Rs particle radius of the LiCoOz electrode (cm) 

S geometric area of the electrode (cm?) 

t time (s) 

19 transference number of the electrolyte 

T temperature (K) 

Ueq equilibrium potential of the electrode (V) 

v thermodynamic factor of the electrolyte 

W active material loading in the electrode (g) 

x x in Li,CoO2, SOC of the electrode 


Greek letters 
Q@,,Q@_ transfer coefficients 
B symmetry factor 


dp, ds electrode or separator thickness (cm) 

£e, &s volume fraction of the electrolyte or active mate- 
rial in solid phase 

Ke conductivity of the electrolyte (S cm!) 

Os conductivity of the solid phase (S cm7!) 

ge, @s liquid or solid phase potential (V) 

Subscripts or superscripts 

c,d charge or discharge 

eff effective value 

n, p positive or negative electrode 

sep separator 


correspond to the C/66, C/13, C/7, C/2.7, and C/1.3 rates, respec- 
tively. For rate capability test at C/66 rate, a single stage constant 
current protocol was used to charge and discharge the half cell. 
For other rates, a two stage constant current protocol was used. 
The cell was first charged or discharged to the desired voltage 
limit (3.0 or 4.35 V) using the target rate. Then a small current 
(C/66 rate) was applied to continue charge or discharge until 
the desired voltage limit was reached again. The use of the two 


Cone shaped Li foil is stuck to 

current collector and placed just 
above the LiCoO> electrode as a 
reference electrode. 


Stainless steel 
current collectors 


LiCoO> 


Separator Li foil 
Fig. 1. Schematic of LiCoO> half cell setup used in the experiment. The refer- 
ence electrode is used to measure the potential of the LiCoO> electrode. 
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Fig. 2. Experimental galvanostatic charge profiles of the LiCoO%2 half cell 
(C=4mA). 


stage charge/discharge protocol was to ensure the cell to reach 
the same SOC at the beginning of charge or discharge for all 
tests. 

The charge and discharge profiles of the LiCoO? electrode are 
presented in Figs. 2 and 3 where the potential is plotted against x 
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Fig. 3. Experimental galvanostatic discharge profiles of the LiCoOz half cell 
(C=4mA). 
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Fig. 4. The electrochemical behavior of the LiCoO> electrode is different during 
charge from that during discharge. The overpotential increases with x during 
galvanostatic discharge, but does not vary much during galvanostatic charge 
(C=4mA). 


in Li,CoO2. The value of x in Li,CoO2 was calculated by using 
the following formula: 


_. @ 
~ 0.274 x W 


where W was the active material loading of the LiCoO3 disc elec- 
trode and was determined to be 0.0245 g from disc electrode size 
and loading data provided by manufacture. The factor of 0.274 
was the theoretical capacity in Ah to convert 1 g LiCoO2 com- 
pletely to CoO2. The accumulative charge or discharge capacity 
Qin Ah was obtained directly from experimental data. The initial 
state (x in Li,CoO2) of the LiCoO>» electrode was determined 
to be 0.99 for charge and 0.393 for discharge from low rate 
potential profiles. 

The potential profiles of the LiCoO2 electrode measured 
during the first stage charge process are compared with those 
measured during the first stage discharge process in Fig. 4 for 
both the C/66 and C/1.3 rates. As seen in Fig. 4, the potential 
profile measured during charge at the C/66 rate almost over- 
lays with the potential profile measured during discharge at the 
same rate, which indicates that the electrode is very close to 
equilibrium. Unfortunately, this is not true for the C/1.3 rate. 
Several observations can be made by comparing the C/66 and 
C/1.3 rate charge and discharge profiles: First, the overpotential 
(relative to the equilibrium potential determined by averaging 
the C/66 rate potentials measured in both charge and discharge) 
does not change much during galvanostatic charge for both the 
C/66 and C/1.3 rates. But it exhibits heavy dependence on state 
of discharge (SOD) during discharge. Second, the potential of 
the LiCoO) electrode starts to drop at high SOD when the dis- 
charge rate is low. But it happens at low SOC if discharge rate 
is high. Third, the potential plateau around x =0.5 (in Li,CoO2) 
gradually disappears in both charge and discharge profiles when 
the current is increased. Many researchers [7—9] explained the 
potential plateau in that region by a phase transition or equilib- 


a) 


= Xo 


rium between ordered and disordered lithium ions in the CoO2 
framework. The movement of large amounts of lithium ions 
in a short time at high current rates may prevent this phase 
equilibrium from occurring, which would explain the gradual 
disappearance of the potential plateau at x=0.5. And lastly, the 
LiCoOy? electrode does not exhibit good rate capability, because 
it reaches the cutoff voltage early at low cell SOD when current 
rate is relatively high (C/1.3). 


3. Model 


Porous electrode models [4—6] have been heavily used in 
the literature to study the performance of lithium ion battery 
systems. A detailed explanation about porous electrode model 
used here can be found elsewhere [5,6]. The model equations 
are summarized below for the convenience of the readers. 

The mass transport in solid LiCoOy particles is described by 
using Fick’s diffusion law: 


aCs 1 ə 2 dCs 
= D; 2 
ot r? or ( sr or ) (2a) 
OC. 
-DE| =0 (2b) 
ər |o 
n 
-D] =4 (2c) 
or r=R, F 


A mass balance on the lithium ions in the solution phase 
yields: 
0 


OE eC, l-t . 
a Fass = 0 (3) 


ot 
where the specific interfacial area of the electrode is as = 3€5/Rs. 


The potential distributions in solid and solution phases are 
described as follows: 


-= V - (DË V ce) — 


V - (@sVġs)—asj=0 (4) 
2RTK dln fx 
ff 0 + 
v(x Voe — - a= (14 A) vina) 
+asj=0 (5) 


The concentrations and potentials in the solid and solution 
phase are coupled together through the intercalation current den- 
sity j which is calculated through the Butler-Volmer equation: 


; : Qa F ; 
J = 10,p | €Xp RT (ds — be Veq JR¢) 


acF ‘ 
exp RT (bs — de Ueq JRe) (6a) 
where the exchange current density io,p is formulated as: 
10,p = kpc? (Cs,max — Cslr=R) Csi R, (6b) 


The computation schematic of a half cell sandwich is shown 
in Fig. 5. The boundary conditions are thus given as: 
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Fig.5. Schematic of a LiCoO; half cell sandwich, consisting of LiCoO2 cathode, 
separator and lithium foil anode. 
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The model equations (Eqs. (2)—(9)) were discretized in space 
using finite difference and solved using the DASSL subroutine 
in Fortran. 

The properties of the LiPFsg based electrolyte, such as 
diffusion coefficient and conductivity are important model 
parameters. But no experimental data were found in the litera- 
ture for the electrolyte system used in this study (1.0 M LiPF6 in 
the EC/PC/EMC/DEC mixture). Valoen and Reimers [10] mea- 
sured the properties of LiPF¢ as a function of temperature and 
concentration in a slightly different electrolyte system which is 
LiPF, in a PC/EC/DMC mixture. Their results were used here 
in the simulation assuming that LiPF¢ behaves similarly in those 


two electrolytes. The diffusion coefficient of the electrolyte was 
found to be [10]: 


4 
4.43 — 2 


7 0.22 x 10°ce 
T —5 x 10°ce — 229 


(10a) 


log(De) = 


The expression for the conductivity of the electrolyte was 
found to be [10]: 


~10.5+0.074T — 6.96 x 10-97? + 668ce : 
Ke = 
© \ -17.8ceT +0.028¢eT? + 4.94 x 105c2 — 886c2T 
(10b) 


The effective value of the diffusion coefficient and conduc- 
tivity of the electrolyte were used in the simulation to account 
for the actual path length of the species in porous media. They 
were calculated through Bruggeman’s correlation: 


DË = pyeBrue (11a) 
Kef — 4 gBrug (11b) 


The thermodynamic factor which account for the non-ideality 
of the electrolyte was found to be [10]: 


v=(1-?) (TE) 


dIn ce 


= 0.601 — 7.59c2° + 3.1 x 104(2.53 — 0.00527)c}5 (12) 


4. Simulation 
4.1. Simulation on charge profiles of the LiCoO2 electrode 
Fig. 6 presents a comparison of the experimental data with 


the predicted charge curve using model equations presented 
above. Clearly, the model predictions do not agree with the 
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Fig. 6. The pseudo-2D model predicts an incorrect potential profile at the begin- 
ning of charge, especially from x = 1.0 to x = 0.8. Symbols represent experimental 
results. Line represents the simulation result. 


Q. Zhang et al. / Journal of Power Sources 165 (2007) 427-435 431 


Table 1 
. Model parameters used to obtain Fig. 6 
Parameter Value Parameter Value 
LiCoOy> electrode 
T (°C) 15? Rg (cm) 10 x 1074? 
W (g) 0.0245? Brug 1.5° 
dp (cm) 64 x 1074? D; (cm? s7!) 2.8 x 107 10¢ 
Sp (cm?) 1.2672 Rg (Qcm?) 200° 
Ee 0.30* X0.p 0.99° 
Es 0.60? Rg (2cm?) 200° 
os (Scm7!) 0.1? kpc (Acm™? (mol cm?) ™!5) 4.2 x 10-2¢ 
19 0.363° 
Separator 
Ceo 1x 10738 Sep (cm?) 1.267? 
ôs (cm) 25 x 1074? Ee 0.46? 
De Eq. (10a) Ke Eq. (10b) 
a-o (14 TË) Eq. (12) 
t dln ce i 


à Manufacture data or experiment data. 
> From Refs. [1-3]. 
€ Values fit to experiment data. 


experimental data for x > 0.8. The model parameters used in 
the simulation are listed in Table 1. The exchange current den- 
sity io» expression and the flat open circuit potential (OCP) of 
LiCoO? from x + 0.75—1 was found to be the combined cause 
of lack of agreement. Fig. 7 shows that the change of io, with x 
during the simulation where it can be seen that io,p has a small 
value at the beginning of charge because cx, max — Cs|r=r, is close 
to zero, which causes the Butler—Volmer equation to predict a 
large overpotential for the set current density. As ig) gradually 
increases during charge because cs max — Cs|r=R, is becoming 
larger, the Butler—Volmer equation predicts a smaller overpoten- 
tial for the same current density. The exchange current density 
in Eq. (6b) was derived based on the assumption of a close rela- 
tion with the Nernstian form of the OCP [4]. However, the OCP 
of LiCoOz around x% 0.75-1 is highly non-Nernstian because 
it shows weak dependence on the lithium ion concentration in 
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Fig. 7. Change of exchange current density ip) with cell SOC when simulating 
Fig. 6. 


the solid matrix (see low rate profiles in Fig. 4). To address this 
problem, the concentration of empty sites Cs max — Cs|r=r, at 
the surface of solid particle is removed from exchange current 
density expression and the io p expression becomes: 


iop = kpce*cslrer, (13) 


Eq. (13) means that the charge kinetics has a relatively strong 
dependence on the concentration of Lit occupied sites in 
LiCoO, but week (or no) dependence on the concentration of 
empty site, similar to kinetics in desorption process. 

Fig. 8 compares the experimental and simulated charge pro- 
files of the LiCoOz electrode using Eq. (13) for io,p instead of 
Eq. (6b). Model parameters used in the simulation are listed in 
Table 2. Simulation results show good agreement with exper- 
iment ones for all charge rates. This study shows that with 
only constant model parameters, the model could predict sat- 
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Fig. 8. Comparison of experiment (symbols) and simulated (lines) charge pro- 
files at different current rates. The rates from bottom to top: C/66, C/13, C/7, 
C/2.7, and C/1.3 (C=4 mA). 
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Table 2 

Model parameters used to obtain Fig. 8° using Eq. (13) 

Parameter Value Parameter Value 
kpc (A cm~? (mol cm?)~!) 6.0 x 1073» Xo.p 0.99» 


à Refer to Table | for parameters not listed here. 
> Values fit to experiment data. 
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Fig. 9. The best simulation results (lines) obtained when model parameters are 
constants. The rates from top to bottom are C/66, C/13, C/7, C/2.7 and C/1.3 
(C=4mA). 


isfactorily the charge profiles of the LiCoO2 electrode where 
overpotential does not vary much with SOC. The model could 
not predict the disappearance of the voltage plateau around 
x=0.5 at high current rates because it does not include any 
mechanism to account for the phenomena. 


4.2. Simulation on discharge profiles of the LiCoO2 
electrode 


Fig. 9 shows the simulation results when the model param- 
eters are kept constant during discharge using Eq. (6b) for io,p. 
The parameters used in the simulation are listed in Table 3. The 
simulated results fit the experiment profiles well at the begin- 
ning of discharge and when current rates are low. But the model 
failed to predict the increased potential drop with cell SOD at 
high rates. Extensive simulations indicate that it is very diffi- 
cult for the model to predict experimental discharge profiles 
when all model parameters are kept constant except those for 
the electrolyte (Eqs. (10) and (12)). 


Table 3 

Model parameters used to obtain Fig. 9° using Eq. (6b) 

Parameter Value Parameter Value 
D; (cm? s7!) 1.4 x 107106 XOp 0.393? 
kp,a (A cm~? (molem™3)~!5) 0.25? 


è Refer to Table 1 for parameters not listed here. 
> Values fit to experiment data. 


In order to achieve a better fit on discharge profiles, an empiri- 
cal expression was introduced to the kinetic expression and used 
along with porous electrode model. The symmetry factor of the 
Li* insertion reaction is assumed to decrease with SOD based 
on our observation of increased potential drop during galvanos- 
tatic discharge. The Li* insertion reaction is usually considered 
as one electron transfer process. 


Lit+S+e =Li-S (14) 


The symmetry factor represents the fractional charge that pro- 
motes the cathodic reaction [4]. When symmetry factor $ 
decreases with SOD, it physically implies that the reversible 
lithium ion intercalation reaction favors more the anodic reaction 
(Li oxidation) than it does the cathodic reaction (Li* reduction). 
That is, it becomes kinetically less favorable to insert Li* into the 
electrode as discharge goes on, which requires increased driv- 
ing force, kinetic overpotential, to maintain the galvanostatic 
discharge current forced by the external circuit. The empirical 
expression for the symmetry factor 6 used is as follows: 


1 
pos ( 1 + exp(a(b — aS os 
A= B ae 
ae (15c) 


Fig. 10 shows the simulation predicted potential profiles dur- 
ing the discharge of the LiCoO 2 electrode at different rates. 
The solid lines represent the simulated profiles and the symbols 
represent experimental data. Model parameters used to obtain 
Fig. 10 are listed in Table 4. Simulation results show good 
agreement with experimental data for all conditions. There is a 
significantly improvement over the previous case (Fig. 9) when 
the symmetry factor is allowed to decrease with SOD. The model 
well captures the fine details of experiment profiles at the begin- 
ning of discharge. The discharge profile of the LiCoO2 electrode 
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Fig. 10. Comparison of experiment (symbols) and simulation (lines) discharge 
profiles at different current rates. Symmetry factor £ is assumed to decrease with 
cell SOD in galvanostatic discharge. 
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Table 4 
Model parameters used to obtain Fig. 10* 
Parameter Value Parameter Value 
D; (cm? s7!) 1.2 x 10710° X0,p 0.393¢ 
kp,a (A cm~? (molcm@3)~!5) 0.25° 
Parameters used in Eq. (15) 
af bt 
Ci66t 35 0.95 
C/13 28 0.94 
C/7 21 0.92 
C/2.7 13 0.88 
C/1.3 8 0.8 


©d Values fit to experiment data. 
1 C=4mA. 
è Refer to Table 1 for parameters not listed here. 


is a Straight line before x reaches 0.5 at C/66 rate. That straight 
line gradually turns into a curve with increased curvature at high 
rate, which could be an indication of increased solid phase diffu- 
sion limitation in the LiCoO? electrode. The model well predicts 
the shape change of the experimental profiles at the beginning 
of discharge using a relatively small solid phase diffusion coef- 
ficient. It does not predict the disappearance of voltage plateau 
around x = 0.5 at high rates for the same reason discussed earlier. 

In order to use the discharge model in other simulations, the 
parameters a and b in symmetry factor expression need to be 
correlated to current rates. The values of a and b obtained from 
fitting experimental discharge profiles are plotted against cur- 
rent rate in Figs. 11 and 12, respectively. Parameter a is found 
to depend rationally on current rate and parameter b linearly. 
The correlation expressions can be found in Table 5. The model 
now could be used to predict discharge profiles of the LiCoO2 
electrode given a discharge current rate. It should be noted that 
the treatment offered in the work, the dependence of symmetry 
factor on the surface concentration, is empirical. Nonetheless, 
the semi-empirical model should yield sufficient good results 
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o Values shown in Table 4 
— Correlation with current rate 
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Fig. 11. The correlation of parameter a with current rate. 
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— Correlation with current rate 
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Fig. 12. The correlation of parameter b with current rate. 


for electrode analysis provided that the operation conditions 
of the electrode are not varied to the extent that the empirical 
representation becomes substantially inaccurate. 


4.3. Simulations using the models for LiCoO2 electrode 


In previous sections, porous electrode models were used to 
simulate experimental charge and discharge profiles of a LiCoO2 
electrode, which is the main purpose of the work. It is shown 


Table 5 
Correlations of parameters a and b for symmetry factor with current rate (xc) 


Parameter Values with 95% confidence intervals 


a=a,l(X_ +2) 


ai 7.453 + 0.841 

a 0.1969 + 0.0285 
b= bi X¢ + bz 

bı —0.2035 + 0.014 

bz 0.9536 + 0.0054 
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Fig. 13. Model predicted charge and discharge profiles for the LiCoO> electrode 
(C=4mA). 


that simulation results agree well with experimental ones. With 
the availability of accurate models, one can study the electro- 
chemical behaviors of the LiCoO) electrode in details. 

The charge and discharge profiles of the LiCoO? electrode are 
simulated under different current rates and the results are shown 
in Fig. 13. The voltage window for the simulation is 3.0—4.35 V. 
The rate capacities are compared to the experiment values in 
Fig. 14. Simulation shows that the LiCoO? electrode used in the 
study does not exhibit a good rate capability in discharge. The 
model predicts the electrode can only achieve about 75% of the 
total capacity at 1 C rate and 50% at 2C rate. 

The rate capacity data can be converted into energy-power 
data to construct a Ragone plot which is a useful tool in com- 
paring different battery designs. It is very time-consuming to 
construct such a plot experimentally. However, computer simu- 
lations can easily and rapidly generate a significant number of 
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Fig. 14. Comparison of experiment and model predicted charge/discharge 
capacity. 
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Fig. 15. Ragone plot of the LiCoO> electrode. 


data points for a Ragone plot. The Ragone plot for the LiCoO2 
electrode is given in Fig. 15. The specific energy (Wh m~?) and 
power (W m~?) are higher than normal values because the elec- 
trode is charged to a higher SOC (x ~ 0.4) which could be the 
most direct way to increase specific energy and power. But in 
real-world applications, the LiCoOz electrode is rarely charged 
to such a high SOC because of the negative effect of phase tran- 
sition at x% 0.5 on the life span of the electrode. The high end 
of charge voltage is used in the rate capability test to help us 
identify the solid phase diffusion limitation in the LiCoO2 elec- 
trode from the shape change of discharge profiles. Ragone plot 
also shows that the electrode is not suited for applications where 
high rate discharge is needed because the specific energy of the 
electrode deteriorates quickly after the current passes 2 C. 

The calculated salt concentration and solution phase potential 
profiles at 2C rate are shown in Figs. 16 and 17, respec- 
tively. The profiles of the salt concentration show moderate 
concentration polarization. And diffusion coefficient of elec- 
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Fig. 16. Calculated solution phase concentration distribution at 2 C discharge 
rate. 


Q. Zhang et al. / Journal of Power Sources 165 (2007) 427-435 435 


0.00 


-0.01 


-0.02 


Electrode 


Electrolyte 


-0.03 


Solution phase potential, ®, (V) 


-0.04 
0.000 0.002 0.004 0.006 0.008 


Position, x (cm) 


Fig. 17. Calculated solution phase potential distribution at 2 C discharge rate. 
trolyte is around 2.5 x 1076 cm? s7! which is several orders 
of magnitude higher than solid phase diffusion coefficient. 
The solution phase potential drop across the cell remains 
near 35mV throughout the discharge. So neither concentra- 
tion polarization nor solution phase potential drop substantially 
limits the cell performance. The solid phase concentration pro- 
files within particles at the electrode—electrolyte interface are 
shown in Fig. 18 for 2C discharge process. The large con- 
centration gradient inside the particle at the end of discharge 
indicates strong solid phase diffusion limitation existing in the 
electrode. 
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Fig. 18. Calculated solid phase concentration profiles in the half cell at 2C 
discharge rate. 


5. Conclusion 


Rate capability study on a LiCoO2 electrode shows that the 
discharge behaviors of the LiCoO) electrode are quite different 
from its charge behaviors. The overpotential of the electrode 
in galvanostatic discharge depends strongly on cell SOD, but it 
does not vary much in galvanostatic charge. 

A pseudo-2D porous electrode model is first used to simulate 
the charge profiles of the LiCoO> electrode. It is found that the 
model predicts false potential profiles at the beginning of charge, 
which is corrected by assuming charge kinetics does not depend 
on the surface concentration of empty sites in LiCoO? particle, 
similar as desorption kinetics. With only constant model param- 
eters, the model predicts potential profiles in good agreement 
with experimental ones. 

Extensive simulations show that the pseudo-2D model has 
difficulty to predict experimental discharge profiles of the 
LiCoO2 electrode when model parameters are kept constant. 
Therefore, symmetry factor in Butler-Volmer equation is empir- 
ically assumed to decrease with cell SOD to account for 
increased potential drop during galvanostatic discharge. It 
means that the insertion of Lit into LiCoO2 electrode becomes 
kinetically less favorable during discharge, which requires an 
increased overpotential to maintain the constant discharge cur- 
rent forced by the external circuit. Simulation results from the 
semi-empirical discharge model are significantly improved and 
match to experimental discharge profiles well. The model should 
yield sufficient good results for electrode analysis given that 
operation conditions of the electrode are not varied to the extent 
that the empirical representation becomes substantially inaccu- 
rate. 
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